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Abstract 

We present first-principles calculations of electronic structures of a class of two-dimensional 
(2D) honeycomb structures of group-V binary compounds. Our results show these new 2D ma¬ 
terials are stable semiconductors with direct or indirect band gaps. The band gap can be tuned 
by applying lattice strain. During their stretchable regime, they all exhibit metal-indirect gap 
semiconductor-direct gap semiconductor-topological insulator (TI) transitions with increasing 
strain from negative (compressive) to positive (tensile) values. The topological phase transi¬ 
tion results from the band inversion at T point due to lattice strain and is irrelevant to spin-orbit 
coupling (SOC). 

Two-dimensional (2D) topological insulators (TIs), also known as quantum spin Hall (QSH) 
insulators, have attracted much attention recent years due to their rich physics and promising ap¬ 
plications in spintronics and quantum computations.- - The QSH effect, which was first proposed 
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in graphene, 5 describes the existence of the edge states on the sides of a 2D TI system. These edge 
states are supposed to consist of two counter-propagating oppositely spin polarized edge channels 
in the band gap of the material. However, the spin-orbit coupling (SOC) in graphene is too weak to 
open a gap large enough to support accessible QSH effect experimentally. Unlike graphene, its el¬ 
emental analogues - silicene and germanene, show energy gaps because they have larger SOC due 
to their buckled honeycomb structures. 6 Therefore, QSH effect can be observed in an experimen¬ 
tally accessible temperature regime in both systems. 7 Other monolayer honeycomb structures of 
group-IV elements as well as III-V binary compounds have been systematically investigated based 
on first-principles calculations-. Theoretical studies show that many of these materials are 2D TIs 
at the equilibrium structure, including stancne-, InBi, GaBi, and TIBi. 1011 Most recently, there 
has been rising interest in layered compounds of group-V elements— -16 . Elemental phosphorus 
have a large number of allotypes. Among them, blue phosphorous has the same layered structure 
as bulk arsenic, antimony, and bismuth. Their corresponding monolayer materials have the same 
buckled honeycomb structure as silicene and germanene. 1217 Monolayer bismuth is reported a 2D 
TI, while other monolayer buckled honeycomb structures of group-V elements are conventional 
semiconductors. 18-20 Furthermore, recent studies indicate that band topology in these 2D materi¬ 
als could be altered by chemical adsorptions, 21-23 external electric field,— ,25 and lattice strain. 25-27 
Under lattice strain, monolayer honeycomb structure of P, As, and Sb all exhibit trivial semicon¬ 
ductor to TI transition. 27 However, the mechanism of the topological phase transition under lattice 
strain has not been revealed. We also note that there is no report of monolayer honeycomb struc¬ 
tures of group-V binary compounds, although bulk AsSb has the same layered structure as bulk 
arsenic, antimony, and bismuth. 28 

In this paper, we study monolayer honeycomb structures of group-V binary compounds (ex¬ 
cept nitrogen) by first-principles calculations. We investigate the stability of these new 2D ma¬ 
terials by studying phonon dispersion and molecular dynamical simulations, and calculate their 
band structures under lattice strains. We find all of them are stable semiconductors with indirect 
or direct energy gaps. Under strain, they all exhibit metal-indirect gap semiconductor-direct gap 
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semiconductor-topological insulator (TI) transitions. It is noteworthy that the mechanism of topo¬ 
logical transition is due to the different response of energies of the valence band and the conduct 
band at T point under lattice strain, and is irrelevant to SOC, which is previously supposed to be 
indispensable for realizing TIs. 

Our first-principles calculations are carried out in the framework of density-functional theory 
(DFT)— within the Perdew-Burke-Ernzerhof generalized gradient approximation (GGA)— imple¬ 
mented in the ABINIT codes. 31 Norm-conserving HGH pseudopotentials 32 and the plane-wave 
cutoff energy of 30 Hartree are used. For monolayer binary compounds, a supercell with a vac¬ 
uum space of 20 A along the z-direction is employed with a 12x12x1 k-point mesh. Both the 
lattice parameters and the positions of all atoms are relaxed until the force is less than 1 meV/A. 
The phonon frequencies are computed using density-functional perturbation theory (DFPT)— with 
a 6 x 6 x 1 q-point mesh. We perform canonical molecular dynamics (MD) simulation at 300K 
with a supercell of 72 atoms. The length of time-step of 3 fs and simulations with 500 steps are 
executed. 

Bulk AsSb, or arsenic antimonide, has the same layered crystal structure as arsenic, antimony, 
and bismuth (space group R3m, No. 166), with the intermediate values of lattice parameters of 
arsenic and antimony, a = 4.025 A, c — 10.84 A. Like monolayer As and Sb, we consider mono- 
layer AsSb consisting of a layer of As and a layer of Sb. This new material has the same structure 
as silicene and germinene, i.e. the buckled honeycomb structure, as shown in Figure |T] We per¬ 
form geometry optimization for such a monolayer AsSb, obtain the lattice constants a = 3.867 
A, and calculate the binding energy. Motivated by this, we consider all possible group-V binary 
compounds (except nitrogen) with the same structure. The optimized lattice constants, buckled 
parameters, and binding energies are calculated, and the corresponding parameters of monolayer 
P, As ,Sb, and Bi are also given for the sake of comparison (see Supporting Information for more 
detail). 

To check the stability, we calculate the phonon dispersion of these monolayer binary com¬ 
pounds. Absence of negative frequencies of the calculated phonon dispersion indicates the stability 
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of monolayer PAs, PSb, PBi, and AsSb (see Figure [2]). Note that although ZA branch (out of plane 
acoustical modes) becomes soft and get imaginary frequencies near T point for AsBi and SbBi, it 
is believed such instability can be removed by the defects. 6-8 On the other hand, the ZA imaginary 
frequencies around T point also depends on the mesh size used in the calculations. It may be an 
artifact of the mesh size. 6-8 The MD simulations also show the stability of these 2D compounds be¬ 
cause no structural collapse happens during the simulations (see Supporting Information for more 
detail). 

Figure [3] shows band structures of these binary compounds. The band structures of monolayer 
phosphorous, arsenic, antimony, and bismuth are also given for comparison (see Supporting Infor¬ 
mation for more detail). The calculated band structures (without SOC) show similarities between 
different systems. This is not surprising, since these elements are chemically similar. Taking into 
account SOC, the number of bands of binary compounds is doubled due to the nonsymmetric struc¬ 
ture and SOC. Note that SOC just modifies the band gaps but does not change band orders for these 
binary compounds, P, As, and Sb. However, the band structure of Bi differs from others when SOC 
is included. The parity of the states that form the valence and conduction bands is reversed at T 
point. This exchange of bands is the consequence of the increased SOC in Bi, and Bi is reported 
a 2D topological insulator 18 . On the contrary, other 2D group-V elemental materials and binary 
compounds with the same buckled honeycomb structure are conventional insulators. Take mono- 
layer AsSb for example, the calculated band structure (without SOC) shows it is a semiconductor 
with a indirect gap of 1.48 eV, while a gap of 1.73 eV at T. The valence band maximum (VBM) is 
at r point, and the conduct band minimum (CBM) is at k = (0.132,0) point. The band structures 
(without SOC) of binary compounds containing Bi possess direct gaps, while others show indirect 
gaps. 

We perform strain engineering to these 2D binary compounds. Figure 0] shows calculated 
energies, pressures, buckling parameters, and energy gaps as functions of the strain. From the place 
where the stress maximum(minimum) occur we can determine the stretchability of the material. 
Take AsSb for example, the stretchability is about —12% for press, and 18% for stretch. During 
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these domain, buckling parameters h and 0 change almost linearly. When strain is larger than 
26%, h drops to zero and 0 equals 90°, and corresponding energy also changes rapidly. It means 
the system becomes planar honeycomb structure as graphene. Band gaps shown in Figure |4(e)| 


(without SOC) and Figure |4(f)| (with SOC) indicate that the system is metallic when strain is 
between —12% and —8%; when strain is between —8% and 2%, the system is a indirect insulator; 
the system becomes a direct semiconductor with a gap at T when strain is larger than 2%. The 
band gap decreases with increasing strain and closes at about 12% strain, then reopen when strain 
is larger than 12%. Band structures of other binary compounds as well as P, As, and Sb under 
strain have similar characteristics (see Supporting Information for more detail). 

To get a better understanding of the mechanism underlying the band inversion, we investigate 
the band structure evolution at T point. Take AsSb for example, As and Sb have five electrons in 
the outer shell. For each atom, two electrons occupy s orbital and three electrons occupy p orbitals. 
They form five energy bands below the Fermi energy. The two s bands are much lower than the 
three p bands. The three p bands of two atoms in the cell form six bands with three below the Fermi 
energy and three above the Fermi energy. According to our first-principles calculations, ignoring 
normalization factors, these six bands around F point are mainly denoted by |1) = \pf) — /;?), 
|2) — | pf) + | p B ) (in fact, there are some (v 4 ) and \s B ) components in |2)), |3) = p 4 ) — |pf), 
|4) = | Py) — \p B ), |5) = |p 4 ) + \p B ), and |6) = |p 4 ) + | p B ). A denotes As atom and B Sb atom. 
The 11), 13), and |4) are below the Fermi energy while 12), |5), and |6) are above the Fermi energy. 
When SOC is not included, |3) and |4) are degenerated, and so are |5) and |6). 

Under ambient conditions, 11) is below |3) and |4). When strain reaches up to 2%, the indirect 
gap becomes a direct gap. When strain is 3%, |1), 13), and |4) have the same energy. When 
strain is beyond 3%, the two bands below and above the Fermi energy at T point is constructed 
by 11) and 12), respectively. Figure [5] shows the band structures of monolayer AsSb at strain of 
8% and 16% respectively, and the band composition. We see the band inversion at T according 
to Figure |5(a)| and Figure |5(c)[ Because h is relatively large when strain is small, the low half of 
|p 4 ) is close to the upper half of | p B ), their coupling is bonding when they are out of phase, while 


5 









their coupling is anti-bonding when they are in phase. Therefore |1) has a lower energy than |2). 
Figure [5(b)] shows 11) and 12), corresponding to VBM and CBM, respectively. In other words, 11) 
is the highest occupied molecular orbital (HOMO) at T point, while |2) is the lowest unoccupied 
molecular orbital (LUMO) at T point. Note that |2) shows [v 4 ) and \s B ) components. With further 
increase of strain, h decreases and the low half of \p£) is move from upper half of \pf) to low half 
of | p?), then energy of 11) increases and |2) decrease. When strain is about 12%, |1) and |2) have 
the same energy and the gap closes. When strain is larger than 12%, the low half of | pf) is closer 
to the low half of \p B ) because of small h, then their coupling is bonding when they are in phase, 
while their coupling is anti-bonding when they are out of phase, therefore 11) has a higher energy 
than 12). Now |2) is HOMO at T point, while |1) is LUMO at T point, and the band inverts, as 
shown in Figure [5(c)] and Figure [5(d)] 

When A and B are the same element such as P, As, and Sb monolayer, the system has inversion 
symmetry. |1) and |2) orbitals have even and odd parity, respectively. The 7Li invariant can be 
determined from the parities of the occupied states at the time-reversal invariant momenta in the 
Brillouin zone. 34 There are four time-reversal invariant momenta for buckled honeycomb structure, 
i.e., one F and three M points, as shown in Figure [1(c)] The band inverts at F means topological 
invariants 7Li changes. Directly parity counting at time reversal invariant momenta shows that the 
Z 2 changes from 0 to 1 for monolayer P, As, and Sb, when stain is larger than 12% (see Supporting 
Information for more detail). We note that As and Sb are reported to exhibit topological phase 
transition at tensile strain 6%.- 27 However, this critical value is estimated according to bulk As and 
Sb lattice constants a , which is about 6% longer than that of corresponding 2D lattice constants. 26 
Taking monolayer As and Sb lattice constants as bases, the critical value should be 12%. 

When A and B are different elements, the inversion symmetry is broken. Then |1), 12), and 
other states, are not the eigenstates of parity. Then it is not possible to determine the 7Li invariant 
from the symmetry of the occupied states at time-reversal invariant momenta. To check the topo¬ 
logical invariants, we adopt the method proposed by Soluyanov and Vanderbilt. — The method is to 
track the largest gap between Wannier charge centers (WCCs). The path following the largest gap 
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between WCCs cross the WCCs bands a number of times that is equal, mod 2, to the 7Li invariant. 
The calculated TLi indices show AsSb is a conventional insulator under 8 % strain and a topological 
insulator under 16% (see Supporting Information for more detail). 

The mechanism of the band inversion at T point is the different response of energies of 11) 
and 12) to lattice strain. The origin of the nontrivial topology in buckled honeycomb structures of 
group-V binary compounds (as well as monolayer P, As, and Sb with the same structure) results 
from the band inversion due to lattice strain. It is neither like the band inversion at T point in 
monolayer Bi, HgTe quantum well, and BhSe 3 , which is due to SOC, nor like that in silicene and 
germanene, which results from massive Dirac cone and there is no band inversion. In the case of 
these 2D honeycomb structures of group-V binary compounds, SOC just modifies energy gaps, 
as shown in Figure |4(e)| and Figure |4(f)[ For example, monolayer AsSb, PBi, and SbBi exhibit 
topological energy gaps of 77, 333, and 69 meV under 16% lattice strain without including SOC, 
respectively. Taking into account SOC, the corresponding energy gaps are 50, 29, and 247 meV, 
respectively. 

Note that the existence of ,v 4 ) and | s B ) components in |2) does not change the mechanism of 
the band inversion mentioned above. However, it makes the difference between band dispersions 
of 11) and 12) near T point. A remarkable feature can be seen that the 11) band shows more flatter 
with increasing strain, in particular when it gets close to the Fermi energy. In fact, the |1) band 
around T point can be approximately given by (see Supporting Information for more detail) 

£(k) = -3 [{V ppa ~ V ppK ) cos 2 e + V ppn \ (M - Bk 2 ) =E\(M — Bk 2 ) (1) 

When the energy of 11) equals to that of |2), we have E\ — 0. It means the valency band becomes 
flat. On the contrary, the existence of \ A ) and \s B ) components in \2) makes its band dispersive. In 
recent years, lattice models with flat bands have attracted attention for a number of reasons, among 
them are enhanced interaction effects . 36 

The existence of topological protected gapless edge states is one of the most important con- 
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sequences of 2D topological insulators. To verify the edge states, we construct zigzag ribbons 
with 48 atoms in a unit cell (about 84.40 A and 90.65 A wide), under strain 8% and 16%, respec¬ 
tively. Figure [6] shows the calculated band structures. We observe both ribbons show band gaps, 
and inside the gap, as expected, we observe four states, corresponding to two spin-splitting edge 
states from either side of the ribbon inside bulk gap. AsSb ribbon has topologically protected gap¬ 
less edge states under 16%. Although gapless edge states are also shown under 8%, they are not 
topologically protected. 

In conclusion, we predict that monolayer honeycomb structures of group-V binary compounds 
are stable according to phonon dispersion and molecular dynamical simulations. These new 2D 
materials are semiconductors under ambient conditions. Monolayer PBi, AsBi, and SbBi possess 
direct gaps, while PAs, PSB, and AsSb possess indirect gaps. Under lattice strain, they show 
topological phase transitions according to calculated Z 2 invariants. By constructing ribbons, we 
show topologically protected edge states exist at nontrivial topology phase. Analysis of the band 
structure evolution at T point indicates that the band inversion is due to lattice strain and irrelevant 
to SOC. This mechanism of topological transition under lattice strain is different from that of 
previously studied 2D TIs, in which SOC is indispensable for nontrivial band topology. Therefore 
it is possible to seek nontrivial band topology in materials composed of constituents with lighter 
atoms and thus smaller SOC. We hope that our work will promote the research aiming at the 
synthesis of these new 2D materials and the search for new 2D TIs. 
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Figure 1: The lattice geometry of buckled honeycomb structure, a: Top view of buckled hon¬ 
eycomb structure. Note that A sublattice and B sublattice (denoted by green and brown color, 
respectively) are not coplanar. a\ = a(l,0), a 2 = a(— A, ^). b: Side view of buckled honeycomb 
structure. Definition of the angle 0 and buckling parameter h. c: Brillouin zone and specific 
symmetry points. b\ = ^(1 &),b 2 = ^(0,^) 
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(a) PAs 


(b) PSb 


(c) PBi 





(d) AsSb 


(e) AsBi 


(f) SbBi 


Figure 2: Phonon dispersion curves of buckled honeycomb structures of group-V binary com¬ 
pounds. 





(a) PAs 


(b) PSb 


(c) PBi 





(d) AsSb 


(e) AsBi 


(f) SbBi 


Figure 3: Energy bands of binary compounds with buckled honeycomb structure. 
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Figure 4: Some parameters as functions of lattice strain, a: Total energy, b: In-plane pressure, c: 
Buckling parameter h. d: Angle 9. e: Energy gap (without SOC). f: Energy gap (with SOC). 



Figure 5: Band structures of AsSb at different lattice strains, HOMO and LUMO orbitals at 
r point, a: Band structure at 8% lattice strain, b: HOMO (|1) = | p^) — \p?)) and LUMO (|2) = 
\Pz) + I Pz) ~ ■ s ' / ') + |s s )) orbitals at F point at 8% lattice strain. Green and brown spheres represent 
As atom and Sb atom, respectively. Yellow (green) color in orbitals means positive (negative) 
phase, c: Band structure at 16% lattice strain, b: HOMO (|2)) and LUMO (|1)) orbitals at T point 
at 16% lattice strain. 
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(a) (b) (c) 

Figure 6: Band structures of AsSb nanoribbon under different lattice strains, a: Band structure at 
8 % lattice strain, b: Band structure at 16% lattice strain, c: The same band structure as shown in b. 
Red edge state comes from As atoms at one edge of the ribbon, while blue edge state comes from 
Sb atoms at the other edge of the ribbon. 
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